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Abstract 

The  primary  spanwise  organized  vortex  structure  and  the 
secondary  streamwise  vortex  structure  of  turbulent  mixing  layers  have 
been  investigated.  Flow  visualization  motion  pictures  of  a  constant 
density  mixing  layer  were  used  to  measure  the  properties  of  the  large 
scale  vortices.  It  was  found  that  after  an  initial  transition  region 
mean  properties  of  the  large  scale  vortices  reach  the  expected  linear 
growth  with  downstream  distance  required  by  similarity.  In  the  self¬ 
similar  region,  the  vortex  core  area  and  visual  thickness  increase 

continuously  during  its  life-span. 

A  theoretical  model  of  probability  distribution  function  for  the 

large-scale  vortex  circulation  was  developed.  This  distribution 
is  found  to  be  lognormal  and  to  have  a  standard  deviation,  normalized 
with  the  mean  of  0.28.  From  this  model  the  mean  life-span  of  the 
vortices  could  also  be  obtained  and  was  found  to  be  0.67  times  the 
mean  life-span  position. 

The  streamwise  streak  pattern  observed  by  Konrad  (1976)  and 
Breidenthal  (1978)  in  plan-view  pictures  of  the  mixing  layer  was 
investigated,  using  flow  visualization  and  spanwise  concentration 
measurements.  It  was  confirmed  that  this  pattern  is  the  result  of 
a  secondary  vortex  structure  dominated  by  streamwise,  counter¬ 
rotating  vortices.  A  detailed  description  of  its  spatial  relation 
to  the  primary,  spanwise  vortex  structure  is  presented.  From  time 
average  flow  pictures,  the  onset  position  and  initial  scale  of  the 
secondary  structures  were  determined.  From  concentration  measurements. 


spanwise  variations  in  mean  properties,  resulting  from  the  secondary 
structure,  were  found.  This  also  showed  an  increase  of  the  spanwise 
scale  with  downstream  distance  and  the  existence  of  the  streamwise 
vortices  in  the  fully  developed  turbulent  region.  In  this  region 
the  mean  spacing  is  found  approximately  equal  to  the  vorticity 
thickness . 
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Preface 


The  turbulent  mixing  layer  has  been  the  subject  of  numerous 
investigations.  Of  particular  interest  are  the  results  of  Brown  and 
Roshko  (1971,  1974)  who  discovered  the  presence  at  high  Reynolds 
numbers  of  spanwise  organized  vortical  structures.  These  vortices  are 
convected  downstream  during  life-times  which  are  terminated  when  they 
interact  or  amalgamate  with  neighboring  vortices  to  form  larger  ones. 

A  form  of  interaction,  the  pairing  interaction,  in  which  two  adjacent 
vortices  rotate  around  each  other  until  they  amalgamate,  was  documented 
by  Winant  and  Browand  (1974)  at  low  and  moderate  Reynolds  numbers. 

This  form  of  amalgamation  and  similar  ones  involving  more  than  two 
vortices  have  also  been  observed  by  us  at  high  Reynolds  numbers  (see 
Roshko,  1976  Dimotakis  and  Brown  1976  and  also  Hernan  and  Jimenez, 

1979) .  The  amalgamation  of  more  than  two  vortices  have  also  been 
studied  in  a  forced  mixing  layer  by  Ho  and  Huang  (1980) .  Moore  and 
Saffman  (1975)  suggested  a  different  possible  form  of  amalgamation 
in  which  a  vortex  is  torn  apart  by  its  neighbors  and  assimilated  into 
them.  Tearing  processes  were  reported  by  Dimotakis  and  Brown  (1976) . 

Winant  and  Browand  proposed  that  the  growth  of  the  mixing  layer 
is  primarily  due  to  amalgamation  processes  but  Hernan  and  Jimenez  (1979) 
who  processed  movies  obtained  by  the  present  author  at  much  higher 
Reynolds  numbers  than  those  in  the  experiments  of  Winant  and  Browand, 
found  that  considerable  (possibly  primary)  growth  occurs  in  individual 
vortices  during  their  lifetime  between  coalescences.  This  observation 
is  consistent  with  the  proposed  model  of  Moore  and  Saffman  in  which 


continual  growth  by  turbulent  entrainment  plays  an  essential  role. 
Their  argument  is  based  on  the  need  to  ingest  irrotational  fluid  into 
the  vortices  so  that  the  necessary  decrease  in  vorticity  can  occur 
as  the  mixing  layer  grows  downstream.  At  low  Reynolds  numbers  this 
ingestion  could  occur  by  viscous  diffusion,  as  in  the  flows  calculated 
by  Patnaik,  Sherman  and  Corcos  (1976),  but  at  high  Reynolds  numbers 
it  would  have  to  occur  by  turbulent  mixing  as  suggested  by  Moore  and 
Saf fman. 

These  views  of  the  primary  vortex  structure  and  processes  in  a 
mixing  layer  illustrate  the  need  to  consider  the  effects  of  what  might 
be  called  residual  turbulence  on  those  processes,  i.e.,  the  effects 
of  secondary  and  possibly  high  order  instabilities  (considering  the 
spanwise  organized  vortices  to  be  the  structures  resulting  from  the 
primary  instability). 

Pertinent  to  this  discussion  is  the  transition  discovered  by 
Konrad  (1976)  and  further  explored  by  Breidenthal  (1978),  who  named 
it  the  "mixing  transition."  The  primary,  spanwise  vortices  and  their 
amalgamations  exist  on  both  sides  of  the  mixing  transition,  which  is 
characterized  mainly  by  the  appearancepn  the  high-Reynolds-number  side 
of  the  transition, of  smaller  scale,  three  dimensional  structures  (as 
seen  on  flow  visualizations)  and  by  increased  "mixedness"  (as  exhibited 
by  the  extent  of  chemical  reaction  between  reactants  carried  by  the 
two  streams).  The  development  of  the  smaller  scale  structure  was 
believed,  by  Konrad  and  Breidenthal,  to  be  associated  with  streamwise 
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oriented  vortex  structures  whose  presence  was  inferred  from  the  stream- 
wise  streaks  that  were  visible  on  plan  views  of  the  mixing  layers. 

In  the  present  investigations  our  purpose  was  to  explore  further 
the  properties  of  the  primary  structures  and  of  the  development  and 
effects  of  secondary,  three  dimensional  structures.  The  results  are 
presented  in  two  parts. 

In  Part  I,  the  statistical  properties  of  the  primary,  spanwise 
oriented  structures  are  obtained  by  examining,  frame-by-frame,  shadow¬ 
graph  motion  pictures  of  the  edge  view  of  mixing  layers.  A  theoretical 
model  of  the  distribution  of  vortex  strengths,  following  the  lead  of 
Takaki  and  Kovasznay  (1978),  is  also  developed  and  shown  to  give  a 
better  fit  than  their  model  to  the  experimental  data.  It  also  provides 
insights  into  the  physical  processes  occurring. 

Part  II  is  devoted  mainly  to  an  exploration  of  the  effects  and 
nature  of  the  secondary  structure  associated  with  the  streamwlse 
streaks  observed  by  Konrad  and  by  Breidenthal.  In  particular,  flow 
visualization  in  a  plane  normal  to  the  flow  direction  confirms  their 
conjecture  that  the  streaks  are  associated  with  streamwlse  vortices. 

By  examination  of  successive  frames  in  a  movie  of  that  plane,  a  three 
dimensional,  spatial  picture  of  the  secondary  structure  and  its 
relation  to  the  primary  structure  is  constructed.  It  Is  consistent 
with  results  from  concentration  measurements  that  were  also  made  in 
a  mixing  layer  between  streams  of  helium  and  nitrogen. 
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1.1  Introduction 

In  this  part,  we  describe  our  experimental  and  theoretical 
studies  of  the  properties  of  the  large  vortices  in  a  turbulent  mixing 
layer.  It  had  been  found  by  previous  investigators  that  in  a  well 
developed  turbulent  mixing  layer  these  vortex  structures  have  a  broad 
distribution  of  sizes  at  any  given  downstream  position.  This 
distribution  is,  in  a  sense,  the  primary  manifestation  of  the  turbulent 
character  of  the  well  developed  flow.  Another  manifestation  is  the 
presence  of  smaller  scale,  three  dimensional  structures  and  motions, 
but  these  play  only  an  indirect  role  in  part  I.  They  will  be  discussed 
more  specifically  in  part  II, 

Brown  and  Roshko  (1974)  measured  the  probability  distribution  of 
vortex  spacing  in  a  mixing  layer  between  helium  and  nitrogen.  Winant 
and  Browand  (1974)  obtained  a  similar  probability  distribution  in 
water  at  moderate  Reynolds  number.  On  the  basis  of  these  results 
Brown  and  Roshko  concluded  that  the  mean  spacing  normalized  with  the 
local  vorticity  thickness  is  constant  independent  of  velocity  and 
density  ratio,  l  *  36^.  In  these  investigations  self-similarity  was 
shown  for  the  mean  velocity  and  density  profiles  but  was  not  shown  for 
the  mean  properties  of  the  vortices.  Self-similarity  of  the  vortex 
spacing  is  implied  by  the  self-similarity  of  spectra  obtained  at  various 
downstream  locations  (e.g.,  Kolpin,  1964)  and  of  autocorrelations 
obtained  by  Dimotakis  and  Brown  (1976)  who  found  3,1  <  1/6  <  5.0, 

U) 

Limited  experimental  data  are  available  on  the  evolution  of  the 
vortices  during  their  life-times,  Brown  and  Roshko  (1974)  found  the 
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spacing  between  vortices  to  be  constant  during  life-time.  On  the  other 
hand,  the  vortex  cores,  as  measured  by  Hernan  and  Jimenez  (1979) 
increase  in  size  during  their  life-times.  Furthermore,  they  found  the 
increase  in  size  of  the  vortex  core  during  a  life-time  to  be 
appreciably  larger  than  during  emaigamation,  i.e.,  little  additional 
growth  of  the  cores  is  accompli:''  id  by  amalgamation. 

An  aspect  of  the  vortex  eve*  on  which  there  had  been  no 

information  so  far  is  the  eftect  of  the  conditions  at  formation  of 
a  vortex  on  its  subsequent  evolution.  In  the  fully  developed 
turbulent  region  where  the  vortex  cores  can  be  considered  turbulent 
the  significant  parameters  are  the  characteristic  time  of  the  small 
scale  motions  within  the  vortex  core  and  the  life-time  of  the  vortex. 

If  the  former  is  much  longer  the  initial  conditions  will  determine 
the  evolution  of  the  vortices.  However,  if  the  characteristic  time 
of  the  small  scale  motions  is  much  shorter,  the  vortex  evolution  can 
be  expected  to  reach  an  asymptotic  state  characterized  only  by  the 
relevant  global  parameters,  namely  its  circulation  and  position. 

In  this  investigation  we  have  used  flow  visualization  motion 
pictures  of  a  constant  density  mixing  layer  to  measure  the  properties 
of  the  large  scale  vortices,  to  follow  their  evolution  and  to 
characterize  their  statistical  properties. 

In  addition  we  have  studied  the  statistical  theory  of  vortex 
circulation  for  the  self-similar  mixing  layer.  This  theory  was 
previously  developed  by  Takaki  and  Kovasznay  (1978)  who  found 
reasonable  agreement  with  experimental  results.  Our  formulation  is 
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somewhat  different  in  that,  among  other  things,  the  results  obtained 
for  the  probability  distribution  of  vortex  circulation  are  also  used 
to  calculate  the  normalized  life  span  of  the  vortices. 

In  this  part  we  first  present  the  experimental  results,  in 
Chapter  1,2,  and  then  the  theoretical  discussion,  in  Chapter  1.3, 
followed  by  a  summary  and  discussion. 
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1.2  Flow  Visualization  Results. 

Two  different  shadowgraph  motion  pictures  of  a  constant  density 
mixing  layer  at  different  Reynolds  number  have  been  used  to  measure  the 
visual  properties  of  the  large  scale  vortices.  The  experimental 
technique  is  described  first.  The  results  concerning  the  vortex 
circulation  (spacing)  follow.  They  were  obtained  at  relatively  low 
Reynolds  number,  so  particular  emphasis  is  made  in  the  distinction 
between  the  transition  region  and  the  fully  developed  turbulent  region. 
In  the  last  section  the  measurements  of  the  visual  thickness  and  core 
area  evolution  are  presented.  These  results  are  compared  with  several 
two-dimensional  models  of  the  evolution  of  the  large  scale  vortices. 
Throughout  this  presentation  reference  will  be  made  to  the  results  of 
the  statistical  theory  presented  in  the  next  chapter. 

1.2.1  Apparatus  and  Data  Processing 

Flow  visualization  motion  pictures  were  obtained  in  the  gas  mixing 
layer  facility  described  by  Brown  and  Roshko  (1974).  Nitrogen  and  a 
mixture  of  33.33%  helium  in  argon  were  chosen  as  working  fluids  because 
their  difference  in  optical  refractive  index  allows  shadowgraph  flow 
visualization  of  a  constant  density  mixing  layer  (Konrad,  1976).  Two 
motion  pictures  were  obtained  at  the  same  nominal  flow  velocities  of 
1000  cm/s  and  380  cm/s,  but  at  different  static  pressures,  4  and  8 
atmospheres.  A  sketch  of  the  mixing  layer  and  coordinate  system 
used  is  presented  in  figure  1.1.  The  parameters  characterizing 
the  evolution  of  the  layer  are  the  velocity  ratio  r  *  0.38  and  the 
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Reynolds  number  AU<5  /v„  =  4.6  x  103  and  2.5  x  104  for  the  4  and 

to  N2 

8  atmospheres  movies  respecitvely .  These  values  of  the  Reynolds  number 
were  calculated  using  an  estimate  of  the  vorticity  thickness  at 
x  =  8  cm,  the  test  section  length  used  in  the  experiments. 

The  shadowgraph  optical  system  is  sketched  in  figure  1.2.  A 
mercury  arc  lamp  (.1  mm  rectangular  light  source)  is  positioned  at  the 
focal  point  of  a  spherical  mirror;  the  resulting  collimated  beam, 

20  cm  in  diameter,  is  sent  through  the  test  section  in  the  spanwise 
direction.  The  beam  is  then  focused  on  the  film  by  a  1200  mm  f-6.3 
lens  and  a  75  mm  camera  lens  located  at  the  focal  point  of  the 
preceding  lens.  A  Hycam  model  K20S4E  high  speed  camera  was  used  at  a 
nominal  frame  rate  of  3000  pictures  per  second.  Actual  frame  rates 
were  measured  using  timing  marks  recorded  on  the  film. 

A  typical  frame  from  one  of  the  motion  pictures  is  shown  in 
figure  1.3.  The  dark  features  within  it  mark  regions  with  large  index 
of  refraction  gradients  on  a  plane  perpendicular  to  the  direction  of 
propagation  of  the  beam,  hence  such  regions  mark  an  interface  between 
fresh  fluid  form  the  two  streams  or  between  mixed  fluid  and  one  of  the 
streams.  Adjacent  vortices  are  separated  by  S  shaped  lines  inclined 
in  the  flow  direction  which  are  called  the  "braids"  following  the 
terminology  of  Patnaik  et  al.  (1976).  These  features  are  convected  by 
the  flow,  with  some  disappearing  during  the  amalgamation  process. 
Pairing  (Winant  and  Browand,  1976)  is  the  most  frequent  form  of 
amalgamation  found  in  these  pictures;  it  involves  the  rotation  of  two 
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vortices  around  each  other  until  the  intervening  braid  can  no  longer 
be  identified. 

For  the  purpose  of  this  investigation  a  large  scale  vortex  is 
defined  as  the  region  between  two  successive  braid/x-axis  intersections 
irrespective  of  the  braid  orientation.  For  example,  the  vortices  in 
the  left  middle  of  the  picture  in  figure  1.3  were  considered  two 
different  vortices  even  though  the  process  of  pairing  may  have  started 
already.  Whereas  the  nearly  completed  pairing  in  the  right  middle  was 
considered  to  be  one  new  vortex.  A  consequence  of  this  definition 
is  that  amalgamation  occurs  between  consecutive  frames  when  a  braid 
disappears . 

The  large  scale  vortex  properties  which  were  measured  are  sketched 
in  figure  1.4;  these  are  the  braid/x-axis  intersection,  the  area  Ac 
and  the  visual  thickness  D  of  the  core.  The  measurements  of  the 
braid/x-axis  intersections  were  made  with  a  Kodak  film  analyser  on  the 
film  obtained  at  p^  =  4  atm.  This  parameter  and  the  area  and  visual 
thickness  of  the  core  on  the  film  obtained  at  8  atm.  were  measured  with 
a  HP  9864  A  digitizer  (section  1.2.3).  In  both  cases  reference 
features  (e.g.,  the  channel  width)  on  each  frame  were  utilized  to 
obtain  the  actual  magnifications.  The  area  of  the  cores  was  measured 
with  the  digitizer  by  moving  the  cursor  along  its  periphery.  The 
accuracy  of  the  measurements  was  limited  by  the  finite  thickness  of  the 
features  that  define  the  large  scale  vortices  on  the  film  itself.  They 
are  associated  in  part  with  the  spanwise  nonuniformities  described  in 
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part  II.  Their  thickness  is  approximately  constant  throughout  the 
length  of  the  layer;  when  scaled  to  the  experiment  size  they  give  a 
measurement  error  of  0.3  cm. 


The  measured  locations  of  adjacent  braid/x-axis  intersections 
were  used  to  calculate  the  distance  2  between  them  (figure  1.4)  and 
the  mid-point  between  them.  The  former  characterizes  not  only  the 
spacing  of  the  vortices  but  also  their  circulation  which  is  given  by 
the  approximate  relation 

r  =  au  •  z 

For  this  reason,  2.  is  called  the  circulation  of  the  vortex  throughout 
this  work.  A  constant  factor  AU  needs  to  be  applied  to  obtain  its 
dimensionally  correct  value.  The  mid-point  between  braid/x-axis 
intersections  is  defined  as  the  downstream  position  of  the  vortex. 

Good  agreement  was  found  between  the  position  of  the  vortex  given  by 
this  definition  and  the  downstream  location  of  the  center  of  mass  of 
the  core  area  (2%  discrepancy) .  It  should  be  noted  that  the  parameters 
measured  are  geometrical  properties  of  the  vortices  which  are  related 
to  changes  in  the  refractive  index  of  the  gas  mixture  and  it  is  not 
clear  how  they  are  related  to  the  vorticity  distribution  within  the 
vortices . 


1.2.2  Vortex  Circulation  Similarity 

The  movie  obtained  at  a  static  pressure  of  A  atm.,  Uj  =  1000  cm/s, 
r  =  0.38,  was  used  to  measure  the  braid/x-axis  intersections  in  a 
sequence  of  600  frames.  From  these  data  the  position  and  circulation  of 
the  vortices  was  calculated. 

A  simple  scheme  was  devised  to  follow  the  vortices  through  several 
frames.  Vortex  realizations  in  two  consecutive  frames  were  considered 
to  be  the  same  vortex  if  their  positions  in  the  two  frames  (xj ,  X2 ) 
satisfied  the  inequality  Ujt  >  x2  -  Xj  >  U2t  where  Vx  ,  V2  are 
the  free  stream  velocities  and  t  is  the  elapsed  time  between  frames. 
When  this  technique  was  applied,  a  total  of  668  vortex  time  histories 
were  obtained,  each  one  giving  the  position  and  circulation  of  the 
vortex  as  a  function  of  time.  Each  time  history  was  approximated  by 
a  straight  line  in  a  least  squares  sense;  the  rms  deviation  from  the 
line  was  found  to  be  always  below  1  mm. 

Self-similarity  was  evaluated  by  considering  the  properties  of 
the  vortices  as  they  were  convected  across  specific  downstream 
positions.  Using  the  time  histories,  the  circulation,  its  time 
derivative  and  the  convection  velocity  of  the  vortices  moving  across 
the  eight  downstream  positions  from  x  =  1  to  8  cm  were  obtained. 

There  is  some  redundancy  in  the  measured  data  which  can  be  used 
to  evaluate  the  vortex  identification  technique.  The  total  length 
convected  through  any  downstream  position  is  given  by  n^*  2.  where 
nT  is  the  total  number  of  vortices  and  2  the  local  mean  circulation. 
This  quantity  is  also  given  by  UcTg  where  Uc  is  the  measured 
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convection  velocity  and  rg  is  the  duration  of  the  sequence  of  600 
frames.  Therefore 


This  result  as  well  as  the  actual  values  of  n^  are  plotted  as 
a  function  of  l  in  figure  1.5.  The  values  calculated  from  the  above 
equation  are  12%  higher  than  those  obtained  by  direct  counting, 
indicating  that  the  technique  fails  to  recognize  some  vortices;  the 
error  is  independent  of  downstream  position.  A  second  feature  shown 
in  this  figure  is  the  reduction  of  the  number  of  vortices  with  down¬ 
stream  position  (increasing  t,) as  the  result  of  amalgamation  processes. 
A  loss  of  110  vortices,  from  165  at  x  =  2  cm  to  55  at  x  -  8  cm, 
occurs  in  the  duration  of  the  sequence. 

The  downstream  evolution  of  the  mean  circulation  is  presented  in 
figure  1.6.  A  linear  dependence  is  found  downstream  of  x^  *  4  cm. 

A  least  squares  fit  to  the  data  downstream  of  this  point  gives 

l  *  .24  (x  -  x0) 

with  xq  «  -  1.4  cm.  This  line  is  also  plotted  in  the  figure.  The 
extent  of  the  transition  region  and  the  location  of  the  virtual  origin 
depend  on  the  initial  conditions,  namely  the  high  speed  boundary  layer 
momentum  thickness  0.  and  Reynolds  number  RQ  (Breidenthal  1978). 
The  boundary  layer  momentum  thickness  for  the  conditions  in  this 
experiment,  calculated  by  Thwaites’  method,  gave 
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These  values  are  compared  in  table  1.1  with  the  results  of  Bradshaw 
(1966)  and  Hussain  and  Zedan  (1978)  in  axisymmetric  mixing  layers  and 
Breidenthal  (1978)  for  a  plane  layer.  The  values  obtained  in  this 
experiment  are  within  the  general  scatter  of  those  results. 

In  the  self-similar  region,  the  growth  rate  of  mean  circulation 
is  dfc/dx  =  0.24.  Using  Brown  and  Roshko's  (1974)  best-fit  line  for 
the  growth  rate  of  the  vorticity  thickness,  d6^/dx  -  0.081,  gives  a 

relation  between  i  and  6  ,  that  is 

U) 

i  =  J.O 

6u> 

This  value  is  in  excellent  agreement  with  the  measurements  of  Brown 
and  Roshko  (1974)  for  a  density  ratio  of  7  and  Winant  and  Browand  (1974) 
for  a  constant  density  mixing  layer,  A  somewhat  larger  value  was 
measured  by  Dimotakis  and  Brown  (1976)  and  by  Koochesfahani  et  al.  (1979). 

The  rms  deviation  of  the  circulation  normalized  with  its  mean  is 
plotted  as  a  function  of  downstream  distance  in  figure  1,7,  A  peak  is 
found  at  x  -  3  cm  in  the  transition  region.  Within  the  self-similar 
region  fluctuations  are  also  encountered.  Their  amplitude  stays  within 
the  estimated  statistical  uncertainty  which,  as  a  result  of  the  small 
sample  available,  is  quite  large  (10%),  The  average  value  of  the 
normalized  rms  deviation  in  the  self-similar  region  (x  >  4  cm)  is  0.24, 
slightly  lower  than  the  values  reported  by  Brown  and  Roshko  (1974)  and 
Winant  and  Browand  (1974).  Koochesfahani  et  al.(1979)  found  a  much 
smaller  deviation,  V/i.  ■  0,03,  due  perhaps  to  their  vortex 
identification  technique. 
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In  the  self-similar  region,  x  >  4  cm,  all  vortices  in  the  sample 
can  be  used  to  obtain  the  normalized  circulation  probability  distribu¬ 
tion.  A  histogram  was  constructed  using  9./ 9.  as  the  similarity 

variable  in  order  to  facilitate  comparison  with  other  experimental  and 
theoretical  results.  This  histogram  is  presented  in  figure  1.8.  Good 
agreement  is  found  with  the  measurements  of  Brown  and  Roshko  and  Winant 
and  Browand.  Also  plotted  in  figure  1.8  is  the  theoretical  result  of 
Takaki  and  Kovasznay  (1978)  and  a  lognormal  distribution  with 
V/l  =  0.28,  the  result  of  the  theory  presented  in  1.3. 

The  vortex  time  histories  (position  x  as  a  function  of  time 
and  circulation  $,  as  function  of  time)  were  also  used  to  obtain  the 
convection  velocity  and  rate  of  change  of  the  circulation  during  the 
life-time  of  a  vortex.  As  we  have  indicated  earlier,  the  variation  of 
position  with  time  was  very  well  approximated  by  a  straight  line  which 
implies  a  constant  convection  velocity  during  each  lifetime.  The  mean 
convection  velocity  of  the  vortices  in  the  selfsimilar  region  is 
724  cm/s  with  standard  deviation  53  cm/s.  The  mean  value  is  5%  higher 
than  the  mean  of  the  free  stream  velocities  690  cm/s.  The  relatively 
large  standard  deviation  results  from  the  error  introduced  by  the 
finite  thickness  of  the  braids.  This  error  can  be  estimated  as 
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where  AUc  is  the  error  in  the  convection  velocity.  Ax  the  error 
in  the  position  (3  mm)  and  x/Uc  the  lifetime  of  the  vortices.  For 
x  *  4  cm,  AUC  becomes  53  cm/s. 

The  rate  of  change  of  the  circulation  during  each  lifetime  in 
the  selfsimilar  region  has  a  mean  value  of  2.5  cm/s  and  standard 
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deviation  of  93  cm/s.  This  surprisingly  large  standard  deviation  may 
be  entirely  associated  with  the  error  from  braid  thickness  which  gives 
an  estimated  error  in  the  circulation  twice  as  large  as  in  the 
position  (i.e.,  106  cm/s).  This  measurement  and  calculation  were  made 
in  order  to  determine  whether  there  are  changes  in  circulation  of  the 
vortices  during  their  lifetimes.  Our  conclusion  is  that,  within  the 
limitations  of  the  measurement,  there  is  no  significant  change  in 
the  mean  value,  but  it  is  not  clear  how  large  any  real  variations  may 
be  because  the  relatively  large  error  from  finite  braid  thickness  tends 
to  obscure  the  results. 

1.2.3  Evolution  During  Life-Span 

A  second  motion  picture  obtained  at  a  static  pressure  of  8  atm. 
and  velocities  of  1000  cm/s  and  300  cm/s  was  used  to  gain  some 
knowledge  of  the  evolution  of  the  parameters  of  the  core  during  its 
life-span.  A  total  of  20  frames  were  processed  in  which  the  vortex 
position  and  circulation  as  well  as  the  core  area  and  visual  thickness 
were  measured.  An  x-t  diagram  of  the  vortices  involved  is  presented 
in  figure  1.9. 

An  estimate  was  made  of  the  distance  required  to  reach  self¬ 
similarity  from  the  results  presented  in  table  1.1.  The  high-speed- 
side  boundary-layer  thickness  and  Reynolds  number  at  separation  are 

0j  ~  .005  cm  Rg  ~  250 


which  give  a  transition  distance  of  1,5  cm,  therefore  all  155  vortices 
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in  Che  sample  (Fig.  1.9)  have  reached  the  self-similar  region.  They 
correspond  to  20  vortex  time  histories. 

The  growth  of  the  mean  circulation  was  obtained  by  fitting  a 
straight  line  through  the  data  in  a  least  squares  sense, 
l  =  0.30  (x  -0.0012) 

Since  dS  /dx  =  .097  from  Brown  and  Roshko  (1974),  1/6  *  3.1. 

0)  w 

The  standard  deviation  normalized  with  the  mean  is  0.27  .  Both 
results  are  in  good  agreement  with  previous  measurements. 

The  core  parameters  D  and  A  (see  Fig.  1.4)  were  measured  for 
each  of  the  155  vortices  in  figure  1.9.  Since  the  circulation  of 
each  vortex  is  constant  during  its  lifetime,  it  seems  reasonable  to 
assume  that  D  and  A  will  be  proportional,  at  least  as  a  first 
guess,  to  l  and  l2  respectively.  Simple  averaging  over  all  155 
samples  gives  (D/2.)  =  0.69  and  (A/ 12)  =  0.49,  with  standard 
deviations  0.16  and  0.18,  respectively.  These  simple  averages, 
however,  do  not  give  information  about  whether  the  cores  grow 
during  their  lifetime.  To  study  this  we  look  for  a  relation  between 
D/l  and  (x  -  Xq  ) / £.  =  £.  Assuming  a  linear  dependence  of  D/2,  on 
£  give  the  best  fit  line 

D/l  =  0.11  q  +  0.31 

with  a  standard  deviation  of  0.08.  This  line  is  shown  in  figure  I. 10. 
Only  three  of  the  20  vortex  time  histories  are  shown,  to  give  an  idea 
of  how  well  the  line  represents  individual  histories.  If  the 
vortices  did  not  grow  during  their  lifetimes,  the  line  would  be 
horizontal  at  the  value  D/l  =0.69.  An  important  point  is  that  the 
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line  of  course  does  not  extend  to  £  =  0  and  £  =  °°,  because  the 
vortices  are  formed  and  disappear  at  finite  values  of  (x  -  x0)/£ 

(see  figure  1.9).  Our  sample  is  too  small  to  obtain  meaningful  mean 
values  for  these  dimensionless  positions  which  we  call  the  mean  forma¬ 
tion  value  and  the  mean  annihilation  value  The  values 

which  are  indicated  by  ticks  on  the  £-axis,  namely  S  «  2.4  and 
*  4.8,  are  deduced  from  the  theoretical  model  discussed  in 
chapter  1.3. 

Instead  of  the  straight  line  in  figure  I. 10,  it  would  be  possible 
to  fit  other  curves,  e.g.,  a  square  root  dependence  with  time 
(i.e.,  with  C  )  as  suggested  by  Moore  and  Saffman  (1975)  to  the  data. 
The  square-root  assumption  gives  D/l  =  0.37^*  with  a  standard 

deviation  of  0.09  . 

The  growth  of  the  vortex  core  area  during  lifetime  was  evaluated 
by  assuming  a  linear  dependence  of  Ac/£2  on  C  resulting  in  the 


relation 


0.13  C 


with  rms  deviation  from  this  line  of  0.11. 


This  line  is  plotted  in  figure  I. 11  along  with  the  same  three 
vortex  time  histories  used  for  comparison  in  figure  I. 10.  Of 
significance  is  the  fact  that  this  best  fit  line  goes  through  the 
origin;  this  result  is  in  agreement  with  Moore  and  Saffman's  suggestion 
for  square-root  growth  rate  (but  here  it  implies  extrapolation  to 
t  *  0) .  It  also  Implies  that  there  is  little  entrainment  into  the 


V 
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vortex  core  during  the  amalgamation  process;  If  we  consider  an 
amalgamation  occurring  at  a  downstream  position  x  involving  two 
vortices  of  circulation  l,  and  £„  and  areas  A  ,  A  ,  the 

12  ci  C2 

newly  formed  vortex  will  have  circulation  i,  ■  fej  +  £  and  area 
Ac  given  by 

A  -  0.13  x  i  *  0.13X  (£,+£)**  A  +  A 
c  '  1  2  c  c 

Brown  and  Roshko  (1974)  used  an  estimate  of  the  rate  of  change  of 

the  core  area  to  calculate  the  entrainment  velocities.  Using  our 

measurements  we  find 


dA 
_ c 

dx 


0.085 


This  value  is  in  relatively  good  agreement  with  the  value  0,077 
obtained  by  Brown  (1974)  from  velocity  measurements. 

The  agreement  found  between  the  entrainment  rate  calculated  from 


the  core  area  measurements  and  that  obtained  from  velocity  measurements 
Indicated  that  the  fluid  engulfed  by  the  large  scale  vortices  eventually 
becomes  part  of  the  mixed  fluid  in  the  core.  At  high  Reynolds  number 
this  process  required  the  presence  of  smaller  scale  three-dimensional 


motions.  In  part  II  a  mechanism  is  described  involving  the  three- 
dimensional  interaction  between  spanwise  and  streamwise  vorticity  which 
results  in  the  growth  of  the  large  scale  core. 

The  vortex  visual  thickness  D  can  be  used  for  comparison  with 
the  numerical  results  of  Patnaik  et  al.  (1976)  who  studied  numerically 
the  nonlinear  development  of  the  Kelvin-Helmholtz  instability  waves 
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at  a  Reynolds  number  of  about  100.  Their  numerical  solution  when 
translated  from  their  temporal  solution  to  a  spatially  growing  shear 
layer  with  r  =  0.3  gives  a  visual  growth  rate  defined  by  their 
density  interface  dSv/dx  **0,10  which  is  in  very  good  agreement  with 
measurement  of  0.11.  On  the  other  hand,  the  maximum  value  of  their 
vortex  core  thickness  is  D/£  ■  0.5,  which  is  about  50%  lower  than 
the  values  observed  in  our  measurements.  Of  course  the  two  situations 
are  quite  different;  in  one  the  diffusion  of  vorticity  in  the  vortex 
and  at  its  boundary  is  by  viscous  (molecular)  mechanisms  while  in 
the  other  it  is  by  (small  scale)  turbulence  mechanisms,  e.g.,  the 
streamwise  vortices  described  in  Part  II. 

Comparison  can  also  be  made  with  some  conclusion  about  core 
parameters  obtained  by  Moore  and  Saffman  (1972,  1975)  and  by 
Pierrehumbert  (1980).  They  found  solutions  of  the  steady  Euler 
equations  for  an  isolated  vortex  in  a  strain  field  and  for  a  row  of 
equal  co-rotating  vortices  respectively,  both  with  uniform  vorticity 
cores.  These  solutions  exist  only  for  values  of  A/ l2  <  0.27  and 
A/ l2  <  0.24  ,  respectively,  where  A  is  the  area  of  the  region 
containing  vorticity.  Our  measurements  however  give  (Ac/Hz)  *  0.49. 
The  discrepancy  may  be  associated  with  the  vorticity  distribution 
within  Ac>  which  may  not  be  uniform,  with  the  streamwise 
nonuniformity  in  vortex  size,  non-steady  effects  or  with  three 
dimensional  effects. 
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1.3  Statistical  Theory  of  Vortex  Circulation 

The  statistical  theory  of  vortex  circulation  attempts  to  describe 
the  self-similar  probability  distribution  from  their  evolution 
properties.  The  implications  of  self-similarity  are  first  introduced, 
they  result  in  a  common  formulation  of  the  Eulerian  (constant  x  )  and 
Lagrangian  (constant  £)  descriptions.  The  properties  of  the  probability 
distribution  are  first  studied  in  an  Eulerian  frame  of  reference  where 
a  conservation  equation  can  be  obtained.  The  results  are  then  applied 
in  the  Lagrangian  frame  of  reference.  The  arguments  presented  here 
are  independent  of  velocity  and  denisty  ratio,  therefore  the  results 
apply  to  any  value  of  these  parameters.  It  also  implies  that  no 
information  can  be  obtained  on  the  growth  rate  of  the  mean  circulation 
since  it  depends  upon  those  parameters. 

1.3.1  Formulation 

The  self-similar  vortex  circulation  distribution  will  be  described 
by  the  probability  of  finding  a  vortex  of  circulation  less  than  l 
at  a  position  less  than  x.  Similarity  requires  it  to  be  a  function 
of  £  *  £/x,  thus 

p(xs  <  x;  £g  <  £)  -  P(V  (1.3.1) 

The  amalgamation  process  will  be  characterized  by  two  probability 
distributions,  the  probability  that  an  amalgamation  will  result  in  the 
annihilation  (P^)  or  formation  (P^,) ,  respectively,  of  a  vortex  of 
circulation  less  than  £  at  a  position  less  than  x.  Using  the 
similarity  scaling  they  can  be  written  as 
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P,(x  <  x,  £  <  £)  =  P.(g) 

As—  s  —  A 

PF(xs  <  x,  £g  <_  £)  *  PF(^) 


(1.3.2) 


Associated  with  each  of  these  probability  distribution  functions 
there  is  a  probability  density  defined  in  the  usual  way 


p(£) 

PA<0 

pf<5> 


dP(£) 

dPA(S) 

dPF(0 


(1.3.3) 


These  probability  density  functions  have  two  different  physical 
interpretations.  For  an  observer  located  at  a  fixed  downstream 
position  x,  in  the  Eulerian  frame  of  reference,  they  describe  the 
probabilities  of  finding,  of  annihilation  and  of  formation  of  a  vortex 
of  circulation  (£,  £  +  d£) ,  respectively.  On  the  other  hand,  for 

an  observer  following  vortices  of  circulation  £,  in  the  Lagrangian 
frame  of  reference,  they  describe  the  probabilities  that  the  event 
will  occur  at  a  position  (x,  x  +  dx) , 

1.3.2  Eulerian  Description 

The  number  of  vortices  of  circulation  (£,  £  +  d£)  being 
convected  through  a  downstream  location  x  per  unit  time  is 

n(£,  x)  -  nT(x)p(£)d£  -  nT(x)5p(£)^  (1.3.4) 

where  n^,(x)  is  the  total  number  being  convected  through  it.  The 


3 


*3! 

•if 


function  n^,(x)  can  be  obtained  by  considering  the  total  length 
convected  per  unit  time,  which  by  definition  equals  the  convection 
velocity  Uc  thus 


nT(x)  =  ^ 


(1.3.5) 


£p(Od5 


If  the  total  number  of  vortices  destroyed  and  formed  by 
amalgamation  per  unit  time  at  a  position  x  are  nA(x)  and  np(x) 
respectively,  the  number  destroyed  and  formed  with  circulation 
0 1,  H  +  d£)  are  given  by 


»a(x)5ps(«  f 


V*>?PF«>  T 


The  combined  effect  of  the  amalgamation  process  is  then 


na(£,  x)d«,  -  [nr(x)  PF(£)  -  nA(x)  PA(Dl£^r 


(1.3.6) 


(1.3.7) 


An  equation  can  be  derived  that  relates  p(£),  p  (^)  and 
p^(£)  by  taking  into  consideration  the  conservation  of  vortex 
circulation  during  its  life-span;  that  is,  the  change  in  n(Jt,  x)  with 
x,  keeping  l  constant,  equals  the  change  due  to  amalgamation 


na(£,  x)  or 


t—  n(£,x)d2,  *  n  (l,  x)  dJl 

dx  a 


m 
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1  +  ^  P(a)  =  ^  8(a)  (1.3.12) 

which  can  readily  be  integrated  to  give 
P  (a)  =  g(a)  -  a+  k 


The  constant  of  integration  k  can  be  determined  by  the  normalization 
condition  on  p(£).  Therefore,  the  general  form  of  p(£)  consistent 
with  similarity  is 


v(V 


eg  (In  O 


CO 


oo 


(1.3.13) 


The  current  understanding  of  the  amalgamation  process  does  not 
allow  a  determination  of  the  function  g(a).  On  the  other  hand,  the 
exponential  dependence  of  p(£)  on  g(£n  £)  suggests  that  the 
properties  of  p(£)  will  be  determined  by  the  behavior  of  g(a) 
around  its  maximum  do >  which,  in  first  approximation,  depends 
quadratically  on  (a  -  cto)  and  therefore  gives  a  lognormal  distribu¬ 
tion  for  p(£),  namely 

p(£)  =  - - -  exp  I - —  (in  §  +  i  (1.3.14) 

/Z?  a£  (  2a2  \  £  /  J 

In  this  expression,  £  is  the  mean  value  of  E,  *  £/x,  while  the 

_ 2  l 

variance  of  £  is  £  (exp  a2  -  1)  where  a2  * - — — r  (Aitchison 

g  1  oi0 ) 

and  Brown,  1957) .  Comparison  with  experimental  results  is  presented  in 
figure  1.8  for  a  value  a  -  .276.  The  agreement  justifies  the  use  of 
this  appproximation  through  the  rest  of  the  discussion. 


V 


T. 
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We  turn  our  attention  now  to  the  width  of  distribution 
characterized  by  the  parameter  a  and  its  relation  to  the  amalgama¬ 
tion  process . 

If  p(£)  given  by  (1.3.14)  is  used  in  equation  (1,3.10)  it  is 
found  that 

kAPA(5)  “  kFPF(?)  *  g2/2  ~  y  ■  P(5>  (1.3.15) 

If  use  is  made  of  the  fact  that  PA(£)  and  Pp(£)  are  positive 
definite,  they  may  be  written  in  the  convenient  form 


pa(U  =  it  [  q2/2  ~  ln  ^  p(^)  H<ep2/2-  VO  +  q(g)J 

A  a2 


(1.3.16) 


=  i r  f  tn-^^r  g-n  p(0  h(c/?  -  e°2/2)  +  q(c)  ] 


where  H  is  the  unit  step  function  and  q(£)  is  a  positive-definite 
function. 

The  functions  p(£),  pA(£)  and  Pp(0  are  sketched  in 

figure  1.12  for  q(£)  =  0.  The  basic  feature  described  by  these 

—  0^/2 

equations  is  that  for  £  <  E,  e  vortices  are  destroyed  by  the 

amalgamation  process  (p.(£)  >  0,  PF(£)  =  0).  The  resulting  vortex 


_  a2  f  2 

formed  by  the  amalgamation  has  £  >  £  e  (pa(£)  “  0»  Pp(£)  >  0)  . 

—  c 7 2 1 2 

It  is  this  transfer  of  vortices  through  £  «  £  e  that  determines 

the  probability  distribution  p(£).  Amalgamations  not  involved  in 
this  transfer  are  described  by  the  function  q(£)  which  does  not 
affect  the  probability  distribution. 


»»'^***-  •c/'.'V- 
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The  normalization  condition  on  pA  and  p?  results  in  the 
equations 


k.  =  1(c)  +  K 
A 

kp  -  I(o)  -  1  +  K 


(1.3.17) 


where 


and 


OO 

'I 


K  =  q(?)  d£ 


I  (a)  5 


/2tF 


+  erf  (a) 


(1.3.18) 


kA  and  kp  are  defined  by  equations  (1.3.9). 


Since 


dnT(x) 

— — —  «*  nF(x)  -  uA(x) ,  kA  -  kp  =  1  and  therefore  the  equations  in 
(1.3.17)  are  not  independent.  The  values  of  kA  and  kp  are 
determined  by  the  amalgamation  process.  In  tearing,  three  vortices 
are  destroyed  and  two  formed,  therefore  kA  =  3  and  kp  =  2; 

3 

similarly.  In  pairing,  kA  =  2,  kp  »  1  and  in  tripling  kA  *  y 
and  kp  »  y  . 

An  additional  assumption  is  needed  concerning  K  in  order  to 
solve  equation  (1.3.17)  for  O.  As  a  first  guess  the  equation 
(1.3.17)  was  solved  for  K  *  0  and  the  above  three  forms  of  amalgama¬ 
tion,  tearing,  pairing,  and  tripling.  These  results  are  given  in 
Table  1.2. 

Using  the  value  a  -  0,276  found  for  pairing,  the  standard 
deviation  £  (exp  a2-  1)^,  normalized  by  £,  is  found  to  be  0.28  in 
good  agreement  with  the  measured  values  of  0.24  for  s  -  1  found  in 
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this  investivation  and  0.31  reported  by  Brown  and  Roshko  (1974)  for  s  *  7. 

This  result  implies  that  coalescence  by  pairing  is  the  dominant 
form  of  interaction,  which  agrees  with  our  observations.  The  presence  of 
other  forms  of  amalgamation  as  well  as  q(£)  ^  0  can  modify  slightly 
the  value  of  <J.  Tearing  results  in  a  smaller  standard  deviation  while 
tripling  gives  a  larger  value  compared  with  the  value  for  pairing.  The 
effect  of  q(£)  ^0  is  to  increase  the  standard  deviation  for  a  given 
form  of  amalgamation. 

The  result  of  this  study  is  that  the  self-similar  probability 
distribution  of  large  scale  vortex  circulation  is  a  lognormal  distribu¬ 
tion  with  a  standard  deviation  of  0,28  independent  of  velocity  and 
density  ratio.  Comparison  is  made  in  figure  1.8  with  experimental 
measurements  in  a  constant  density  mixing  layer  presented  in  the 
previous  chapter  and  for  a  density  ratio  s  =  7  made  by  Brown  and 
Roshko  (1974).  The  agreement  is  excellent,  in  particular  the  location 
of  the  peak,  i.e.,  the  mode  of  the  distribution,  is  well  characterized 
by  this  theory.  Also  plotted  in  the  figure  are  the  results  of  Takaki 
and  Kovasznay  (1978),  who  found  a  wider  distribution  with  the  mode  at 
a  smaller  value  of  H/Jl. 

1.3.3  Lagrangian  Description 

The  results  obtained  in  the  previous  section  provide  some  informa¬ 
tion  on  the  Lagrangian  description  of  vortex  evolution.  The  different 
probability  density  functions  will  be  presented  in  terms  of  the 
similarity  variable  t;  *  x/il  used  in  describing  the  vortex  parameters 


in  section  1.2,3. 


‘■At M 


29 


The  probability  of  finding  a  vortex  of  circulation  2.  at  a 
position  (x,  x  +  dx)  is  given  by  p(Q  d£  which,  after  change  of 


variable,  becomes 


Pc<OdC 


exp  -  —  {In  Vt,  +  ^-)2  d£  (1.3.19) 

2tt  o  2a2 


—  a2  —  —  -3a2 /2 

with  Z,  =  e  /£•  Its  maximum  value  is  found  at  =  £e 

H 

The  probability  of  formation  Pp(£)  and  annihilation  PA(C)  after 
change  of  variable  and  neglecting  q (C )  read 

PF<Od?  =  -  3°2-2  +-ln  L>L-  P,(0  d?  0  <  c/«f  <  e-3°Z/2 
F  [1(a)  -  l]a2  5 


PAU)dS  - 


3a2 /2  +  in  c/i 
I(a)a2 


PC(C)  dc 


e-3°2/2<  c/C  < 


(1.3.20) 


p  (S),  pp(£)  and  p.(£)  are  sketched  in  figure  1.13.  The 

interpretation  of  these  probability  distributions  is  readily  found. 

If  we  consider  vortices  with  a  given  circulation  l,  Pp(G)  indicates 
that  they  are  formed  by  coalescence  until  a  downstream  position 
x  =  Z^SL  is  reached.  From  that  point  amalgamation  gradually  destroys 
the  vortices  of  this  circulation  as  indicated  by  PA^C) • 

The  functions  Pp(5)  and  pA(£)  can  used  to  obtain  the  mean 
formation  and  destruction  locations,  £p  and  ?A,  defined  as 


n- 

'  L. 


it •  ■ 


-e  ~3o2  /2 


C  P,p(0  d?  “  2  5 


-  I (2a)  -  1 


Ko)  -  1 


(1.3.21) 


?  pa<0  dS 
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?  Ko ) 


-  -3o2/2 
Ce 


where  the  function  1(a)  is  given  by  equation  (1.3.18).  The  values 
of  Cp/c;  and  £  /£  for  the  three  forms  of  amalgamation  cosidered 
are  given  in  table  1.2,  as  well  as  the  mean  life-span  L  of  the 
vortices  normalized  with  their  mean  life-span  position  X, 


C  .  C 

_A _ F 


(1.3.22) 


The  value  for  pairing  L/X  *  0.67  can  be  compared  with  the 
measurements  of  Brown  and  Roshko  (1974),  L/X  -  0.39,  and  of  Hernan 
and  Jimenez  (1979),  L/X  -  .61.  The  theoretical  result  overestimates 
slightly  the  experimental  values,  the  presence  of  other  forms  of 
amalgamation  and  some  overlap  may  account  for  some  of  the  discrepancy 
but  it  should  also  be  pointed  out  that  the  experimental  samples  were 
small. 

Comparison  can  also  be  made  with  the  experimental  results 
presented  in  figure  I. 10  and  I. 11  of  the  vortex  evolution.  For  that 


ISM****:  mum*m**  *****  *'*•*■•  ^ 


sample,  C  “  0.3  and  a2  *  .07  which  gives  C  “  3.6,  ^  =  2.4 

and  5.  ■  4.8.  These  mean  values  are  consistent  with  the  measured 
A 

evolution  of  the  three  vortices  presented  in  the  figures. 


I. A  Snmmar 


In  this  investigation  we  have  used  the  properties  of  an  ensemble 
of  large  scale  vortices  to  characterize  the  transition  region  and  the 
subsequent  self-similar  behavior  of  the  flow.  A  model  of  the  statisti 
cal  distribution  of  vortex  strength  (circulation)  was  also  developed 
that  provides  additional  insight  into  the  dynamics  of  the  self-similar 
evolution. 

A  summary  of  the  results  follows. 

(1)  The  extent  of  the  transition  as  determined  in  this  investi¬ 
gation  on  the  basis  of  vortex  similarity  is  in  good  agreement  with 
previous  measurements.  In  the  transition  region  the  mean  circulation 
increases  monotonically  while  the  rms  reaches  a  peak  within  the  region 

(2)  In  the  fully  developed  turbulent  region  the  mean  circulation 

normalized  with  the  vorticity  thickness  is  2/6  *  3.0  and  the  rms 

w 

fluctuation  normalized  with  the  mean  is  2'/2  *  0.24. 

(3)  The  evolution  of  the  normalized  core  parameters,  the  visual 
thickness  D/2,  and  area  A Jt2  are  adequately  described  in  terms  of 
the  similarity  variable  (x  -  xQ)/2  =  £.  The  scatter  of  the  measure¬ 
ments  was  too  large  to  choose  between  a  square  root  dependence  of 
D/2  on  £  or  a  linear  one  (with  non-zero  origin);  the  former  gives 

7  -  0.37  c*f  —  -  0.13  C 

%  t2 

This  functional  dependence  was  suggested  by  Moore  and  Saffman  (1975) . 
The  rms  deviations  from  these  lines,  13%  for  (D/2)  and  22%  for 
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(Ac/g.2  )  respectively,  are  about  the  same  as  the  resolution  of  the 
measurements.  That  is,  any  possible  effects  of  initial  conditions 
(at  formation  of  the  vortices)  are  smaller  than  the  resolution  of  the 
measurement. 

(4)  The  statistical  theory  of  vortex  circulation  presented  in 
Chapter  1.3  adequately  describes  the  shape  and  variance  of  the 
probability  distribution.  It  is  a  lognormal  distribution  with  a 
parameter  a  which  is  independent  of  velocity  and  density  ratio.  The 
assumptions  upon  which  this  result  was  obtained  are 

i.  Self-similarity.  The  probability  distributions  are  only  a 
function  of  the  similarity  variable  £  =  £/x. 

ii.  The  vortex  circulation  is  conserved  during  its  evolution, 

iii.  The  dependence  of  the  amalgamation  process  on  the  circulation 
of  the  vortices  involved  is  such  that  the  function  g(a) 
can  be  approximated  by 

g(a)  =  — (a  -  a0)  +  a(a  -  +  b 

iv.  The  overlap  function  q(£)  =0.  It  implies  that  for  a  given 
value  of  £  amalgamation  can  either  form  or  destroy  vortices 
but  not  both. 

Because  of  the  third  assumption  the  lognormality  of  the  distribu¬ 
tion  should  formally  be  interpreted  as  a  first  term  in  an  expansion 
with  small  parameter  e  -  a3g'"(a0). 

The  numerical  value  of  a  depends  on  the  type  of  amalgamation 
assumed;  if  pairing  is  the  dominant  form  then  o  -  0.276,  which  agrees 
best  with  the  experimental  value. 
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(5)  The  statistical  theory  also  gives  the  mean  formation  and 
annihilation  positions  of  the  vortices.  These  results  are  consistent 
with  the  limited  evidence  available.  The  mean  life-span  normalized 
with  the  mean  position  L/X  =0.67  is  in  reasonable  agreement  with 


experimental  measurements. 
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SECONDARY  STREAMWISE  VORTEX  STRUCTURE 
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II . 1  Introduction 

The  measurements  reported  in  this  part  were  motivated  largely  by 
some  results  of  earlier  experiments  made  at  GALCIT,  notably  those  of 
Konrad  (1976)  and  Breidenthal  (1978),  both  of  whom  noted  on  plan  views  . 
of  the  mixing  layer  the  presence  of  streamwise  streaks  which  were 
clearly  associated  with  the  development  of  three-dimensional  motions. 
Both  of  them  also  described  a  mixing  transition  region  in  which  the 
internal  mixedness  of  the  shear  layer  increases  and  which  was 
associated  with  the  development  of  the  streak  pattern.  It  was 
speculated  that  the  streaks  delineated  streamwise  vortices  which 
contributed  to  the  increased  mixing,  either  directly  and/or  after 
further  instability. 

The  experiments  described  in  this  part  were  done  in  th  2  gas 
Brown-Roshko  facility  used  by  Konrad  as  well  as  in  Breidenthal 's 
water  facility.  It  was  our  objective  to  obtain  a  beter  understanding 
of  the  origin  and  nature  of  the  structure  associated  with  the  stream- 
wise  streak  patterns.  We  approached  this  by  first  observing  the 
effects  of  apparatus  and  flow  parameters  on  the  streak  pattern  early 
in  its  history.  Screens  in  the  settling  chamber  of  the  gas  mixing 
facility  and  the  length  of  the  splitter  plate  were  changed;  Reynolds 
number  and  velocity  ratio  were  varied.  Time  exposure  pictures 
indicated  a  spanwise  mean  structure  which  was  then  investigated  further 
by  means  of  concentration  probes.  These  clearly  indicated  that  indeed 
the  streaks  are  associated  with  streamwise  vortices.  Finally,  direct 
visualizations  of  those  vortices  were  obtained  in  the  water  facility 
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and  high  speed  movies  made  it  possible  to  obtain  a  description  of 
their  spatial  relation  to  the  primary  spanwise  vortices.  The 
implications  of  changes  of  scale  observed  in  the  concentration 
measurements,  and  evident  on  the  earlier  pictures  of  Konrad,  are 
discussed . 
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II. 2  Helium-Nitrogen  Mixing  Layer  Results 

The  Brown-Roshko  (1974)  gas  facility  was  used  to  study  the 
three-dimensional  structure  of  the  mixing  layer  between  helium  and 
nitrogen  streams.  Flow  visualization  and  concentration  measurements 
were  the  principal  experimental  techniques  utilized  in  this  investiga¬ 
tion.  A  brief  description  of  the  apparatus  and  experimental 
techniques,  followed  by  the  flow  visualization  results  and  concentration 
measurements,  will  be  presented. 

II. 2.1  Apparatus  and  Experimental  Techniques 

The  gas  mixing  layer  facility  is  a  two-channel  pressurized,  blow¬ 
down  tunnel.  The  gases  are  stored  in  high-pressure  bottles,  eight 
bottles  for  each  channel.  Two  independent  lines,  which  include  a 
pressure  regulator,  a  metering  valve  and  a  remotely  operated  ball 
valve,  provide  an  adjustable  mass  flow  for  each  channel.  The  working 
section  is  enclosed  in  a  high  pressure  vessel,  whose  pressure  is 
maintained  at  the  required  level  by  a  regulator. 

In  each  channel,  the  gas  flow  enters  an  entrance  section  which 
contains  several  turbulent  reducing  elements,  namely  a  0.64  cm  thick 
foam,  a  perforated  plate  with  44%  open  area  ratio,  a  4  cm  long 
honeycomb  section  with  0,2  cm  cell  size  and  a  fine  mesh  screen  with 
50%  open  area  ratio.  After  that  the  flow  enters  a  two  dimensional 
contraction  with  4  to  1  area  ratio.  The  resulting  free  stream 
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turbulence  level  is  below  0.5%.  The  two  channels  are  separated  by  a 
splitter  plate  which  terminates  at  the  entrance  to  the  test  section 
in  a  6  degree  wedge. 

The  test  section  is  rectangular  in  cross-section,  5  cm  high  by 
10  cm  wide.  At  the  entrance  each  free  stream  occupies  one  half  of 
the  cross-section,  i.e.,  10  cm  x  2.5  cm,  the  larger  dimension  being 
spanwise.  The  top  and  bottom  walls,  those  parallel  to  the  mixing 
layer,  are  adjustable  so  that  the  pressure  gradient  in  the  flow 
direction  can  be  minimized.  The  length  of  the  test  section  is  25  cm, 

but  only  the  initial  11  cm  was  used  in  the  measurements. 

Shadowgraph  and  schlieren  flow  visualization  was  used  in  these 

experiments.  The  optical  system  is  sketched  in  figure  II. 1.  It  is 
based  on  the  system  described  in  Part  I.  Simultaneous  side  and  plane 
views  were  obtained  with  two  first  surface  flat  mirrors  as  shown  in 
the  figure.  In  the  schlieren  system  the  knife  edge  was  positioned 
parallel  to  the  flow  direction,  thus  the  intensity  pattern  on  the  film 
is  proportional  to  the  density  gradient  across  the  layer  in  the  side 
view  and  along  the  span  in  the  plan  view,  respectively.  A  stroboscopic 
light  source,  the  Genrad  1538A,  was  used  for  the  schlieren  pictures. 
Time  exposures  were  obtained  by  strobing  the  unit  at  100  Hz  for 
2  seconds . 

Concentration  measurements  were  made  using  the  probe  developed 
by  Brown  and  Rebollo  (1972).  The  probes  were  calibrated  in  a 


40 


different  facility  against  a  mixture  of  known  composition.  Power 
spectrum  measurements  of  the  concentration  fluctuation  showed  no 
appreciable  energy  content  above  15  kHz.  At  the  velocities  used,  this 
value  corresponds  to  a  spatial  resolution  of  0.2  tnm. 

Two  concentration  probes  and  one  pitot-static  probe  deployed  at 
fixed  positions  along  the  span,  were  used  in  most  data  acquisition 
runs.  During  each  run  the  set  of  three  probes  was  traversed  normal 
to  the  shear  layer,  stopping  at  21  points  within  the  layer.  The 
position  accuracy  was  0.02  cm.  At  each  stopping  point  the  probes' 
output  was  digitized  sequentially  at  a  combined  rate  of  96  kHz  for 
0.18  seconds  and  stored  on  digital  tape  for  later  processing.  The 
data  acquisition,  and  digital  operations  were  performed  with  a  computer- 
controlled  data  acquisition  system  called  the  "Solo  System." 

The  mean  and  rms  fluctuations  of  the  concentration  was  obtained 
from  these  data.  In  addition,  molecular  mixing  was  characterized 
by  computing  the  amount  of  product  formed  in  the  layer  by  an 
hypothetical  infinitely  fast  irreversible  chemical  reaction.  Under 
the  assumption  of  equal  diffusion  coefficients  of  all  the  species 
involved  the  amount  of  product,  normalized  by  the  concentration  of  the 
reactant  carried  by  the  low  speed  fluid,  is  (Toor  1962,  Konrad  1976) 
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where  C,  C  ,  C„  and  C,  are  the  helium,  product  and  free  stream 
P  "0 

reactant  concentrations,  with  f  =  /CA  .  A  single  value  of  f  =  3 
was  used  in  this  investigation  to  characterize  molecular  scale  mixing. 
This  value  was  also  used  by  Konrad  (1976),  hence  comparison  can  be 
made  with  his  detailed  measurements  at  a  fixed  spanwise  position. 

II. 2. 2  Flow  Visualization  Results 

Quantitative  data  on  the  onset  of  the  streamwise  vortices  was 
obtained  from  flow  visualization  pictures.  Figure  II. 2  shows  two 
visualizations  of  a  mixing  layer  with  =  496  cm/s  ,  r  =  0.5, 

p  =  4  atm.  (a)  shows  an  instantaneous  (1-2  psec)  spark  schlieren 
picture,  while  (b)  is  obtained  from  a  2  second  exposure.  It  may  be 
seen  that  the  position  of  the  streamwise  streaks  in  the  two  pictures 
is  the  same,  implying  that  there  is  little  fluctuation  in  their 
spanwise  positions. 

The  persistence  of  the  streaks  at  fixed  spanwise  positions 
suggests  that  they  may  be  associated  with  disturbances  originating 
somewhere  upstream  in  the  settling  chamber.  Indeed,  changes  in 
the  screens  produced  changes  in  their  spanwise  positions.  On  the 
other  hand,  it  is  clear  that  the  flow  disturbances  associated  with 
the  streaks  are  being  amplified  from  their  upstream  origins  until 
they  reach  levels  at  which  they  become  visible  on  the  photographs. 

The  location  of  this  visible  threshold  amplitude  could  be  shifted 
upstream  or  downstream  by  raising  or  lowering  the  Reynolds  number 
(for  fixed  velocity  ratio).  These  flow  changes  also  changed  the 
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number  of  screaks.  For  example,  figure  II. 3  shows  pictures  obtained 
under  exactly  the  same  experimental  conditions  as  figure  II. 2,  except 
for  a  change  of  velocities  as  indicated.  It  may  be  seen  that  there 
is  a  considerable  change  in  the  streak  pattern.  Our  conclusion  from 
these  and  other  comparisons,  discussed  below,  is  that  the  streaks 
are  associated  with  true  flow  instabilities  which,  of  course, 
originate  from  the  available  upstream  disturbances,  but  which  respond 
differently  to  different  flow  conditions.  In  particular,  we  studied 
the  dependence  of  the  mean  spacing  and  the  beginning  of  the  visible 
pattern  on  the  flow  parameters. 

An  appropriate  reference  length  for  the  spanwise  spacing  of  the 
streaks  might  be  the  initial  thickness  of  the  shear  layer  after  it 
leaves  the  splitter  plate.  Apart  from  the  difficulty  of  measuring 
this  small  thickness,  there  are  well  known  problems  of  defining  the 
position  of  initial  thickness,  compounded  here  by  the  fact  that  the 
shear  layer  is  generated  from  initial  boundary  layers  in  two  different 
gases.  Therefore,  we  use  for  the  reference  length  the  wave  length  of 
the  initial  Kelvin-Helmholtz  instability  wave  ^Di  can  easily 

be  measured  on  the  photographs.  In  a  sense,  this  instability  wave 
educes  the  appropriate  length  parameter  and  provides  a  convenient 
means  of  measuring  it. 

Using  the  time  exposures,  measurements  were  made  of  the  distance 
x^D  to  the  visual  origin  of  the  pattern,  of  the  spacing  A^  and 
of  the  visual  thickness  of  the  layer  at  x^.  Instantaneous 
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pictures  were  used  to  measure  the  streamwlse  (Kelvin-Helmholtz) 

wavelength  ^Di  t*ie  rolled-up  laminar  shear  layer  and 

also  A2D  at  x^.  (A2D  is  a^ways  greater  than  ^2Di  because 

amalgamations  have  occurred  before  the  onset  of  the  visible  pattern. 

The  values  of  both  lengths  are  actually  average  values  from  two 

realizations  of  the  flow.)  The  results  of  the  measurements  are 

summarized  in  table  II. 1,  where  the  layer  thickness  at  the  beginning 

of  the  pattern  is  characterized  by  what  we  call  the  onset  Reynolds 

number  AU6  /vN  ,  while  the  spanwise  mean  spacing  and  its  rms 
v0  w2 

fluctuation  are  normalized  by  A2Di  and  X2D  . 

These  results  are  presented  in  figures  II. 4  and  II. 5  for  the  onset 
Reynolds  number  and  the  mean  spacing,  respectively.  We  have  used  in 
these  figures  the  parameter  AU/UC  instead  of  the  velocity  ratio  r  , 
where  AU  Is  the  difference  in  free  stream  velocities  and  Uc  is 
the  convection  velocity  of  the  large  vortices.  Its  relation  to  the 
velocity  and  density  ratio  is  according  to  Brown  (1974) 

AU  =  (1  +  3)(1  -  r) 

Uc  1  +  sr 

As  shown  in  figure  II. 4  the  onset  Reynolds  number  depends  strongly 
on  this  parameter  changing  by  a  factor  of  4  over  the  range  0.8  to  2.2, 
which  corresponds  to  a  change  in  r  from  0.15  to  0.5.  The  extension 
of  the  partition  resulted  in  an  increase  of  the  onset  Reynolds 
number  for  the  same  velocity  ratio.  We  believe  that  this  is  due  to 
a  change  in  boundary-layer  profile  shape. 


The  mean  spanwise  spacing  X »  normalized  with  the  streamwise 

wavelength  at  onset,  A-p  ,  is  plotted  as  a  function  of  the  parameter 

uo 

AU/U^  in  figure  II. 5.  The  error  bars  give  the  rms  variation 
calculated  using  all  the  data  at  a  particular  value  of  AU/Uc>  No 
significant  trend  with  AU/Uc  is  apparent  and  therefore  the  mean  and 
rms  were  calculated  using  all  samples  available,  which  gives 
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The  large  rms  deviation  is  not  surprising  because  individual  realiza¬ 
tions  have  rms  variations  as  high  as  0.3  (see  table  II. 1). 

The  time  exposure  pictures  show  significant  downstream  deteriora¬ 
tion  of  the  coherent  streamwise  pattern  at  sufficiently  high  Reynolds 
number,  as  in  figure  II. 3.  In  these  cases  the  evolution  is  character¬ 
ized  by  a  gradual  loss  of  contrast  and  by  spacing  changes.  Additional 
evidence  of  this  was  obtained  from  high  speed  shadowgraph  motion 
pictures.  They  showed  that  the  spanwise  location  of  the  streaks  is 
initially  fairly  well  fixed,  but  further  downstream  variations  in  the 
pattern  develop.  The  lateral  positions  of  the  streaks  at  a  given 
downstream  location  fluctuate  in  time  and  the  number  of  streaks  across 
the  span  decreases  with  increasing  distance  downstream.  This  behavior 
is  consistent  with  the  time  exposure  pictures  since  the  lateral 
displacement  and  spacing  changes  will  result  in  loss  of  the  initial 
coherence  of  the  pattern.  In  all  flow  conditions  considered  in  this 


investigation  the  streamwise  streaks  were  observed  in  instantaneous 
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pictures  even  at  the  downstream  end  of  the  test  section,  downstream 
of  the  mixing  transition,  up  to  the  highest  Reynolds  number. 

II. 2. 3  Concentration  Measurements 

The  streamwise  streak  pattern  evident  in  the  time  exposures 
indicated  the  existence  of  persistent  spanwise  concentration  gradients 
in  the  spanwise  direction  indicating  that  there  are  spanwise  non¬ 
uniformities  in  the  mean  flow.  The  concentration  measurements  were 
undertaken  to  obtain  quantitative  definition  of  these  nonuniformities 
and  to  provide  some  insight  into  the  underlying  structure  that 
generates  them.  Also  we  wished  to  see  whether  any  mean  spanwise 
structure  actually  remains  in  regions  downstream  of  where  it  has 
disappeared  from  the  time  exposure  pictures.  The  measurements  were 
limited  to  the  case  with  free  stream  velocity  of  1000  cm/s  on  the 
helium  side,  380  cm/s  on  the  nitrogen  side  and  test  section  pressure 
of  6  atmospheres.  These  conditions  s  ■  7  and  r  *  0.38  were  chosen 
so  that  comparison  could  be  made  with  the  results  of  Brown  and 
Roshko  (1974)  and  Konrad  (1976). 

In  order  to  establish  a  basis  for  comparison  with  those  previous 
measurements,  traverses  were  made  with  a  single  probe  positioned  at 
midspan.  The  concentration  thickness  defined  as  the  distance  across 
the  layer  between  points  where  the  mean  concentrations  are  0.01  and 
0.9  was  found  to  follow  the  expected  linear  dependence  on  downstream 
distance 

6  .  -  0,22  (x  -  xB),  x0  -  -  1.58  cm 
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The  growth  rate  0.22  is  in  good  agreement  with  the  results  of 
Brown  and  Roshko  (1974).  The  value  of  xQ  is  used  to  define  the 
similarity  variable  n  =  y/ (x  -  x0 ) • 

Once  the  correspondence  with  earlier  results  had  been  established, 
we  proceeded  to  measure  spanwise  variations  by  using  two  concentration 
probes,  one  of  them  positioned  at  midspan  (  z  »  0)  and  the  other  at 
various  spanwise  positions  z.  The  downstream  evolution  of  the  mean 
concentration  profiles  differs,  depending  upon  spanwise  position.  With 
the  two  probes  spaced  0.64  cm,  pairs  of  profiles  were  obtained  at  four 
downstream  stations  as  shown  in  figure  II. 6.  The  spacing  was  chosen 
on  the  basis  of  time  exposure  pictures  so  that  both  probes  were 
located  at  neighboring  streaks,  consequently  the  close  agreement 
found  at  x  =  3  cm  is  to  be  expected.  As  the  probes  were  moved 
downstream,  keeping  their  spacing  constant,  differences  in  the  mean 
concentration  profiles  developed.  As  will  be  seen  later,  this  is  due 
to  an  increase  in  the  mean  spacing  of  the  streaks  and  they  no  longer 
correspond  to  the  probe  spacing. 

To  study  this  more  completely,  the  mean  concentration  and  its 
rms  fluctuation  were  mapped  over  part  of  the  cross-sectional  (y-z) 
plane,  at  two  downstream  locations,  x  =  4  cm  and  x  -  9  cm.  The 
results  are  plotted  in  figures  II. 7  and  II. 8,  respectively.  At  the 
upstream  location,  the  spacing  of  the  peaks  (or  of  the  valleys) 
corresponds  to  the  streak  spacing  in  the  time  exposure  pictures.  For 
the  mean  concentration,  the  peaks  are  in  phase  (spanwise)  all  across 
the  layer  while  for  the  rms  fluctuation  peaks  on  the  helium  side  are 
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in  phase  with  valleys  on  the  nitrogen  side.  For  the  latter  the  change 
over  occurs  at  approximately  n  =  -  0.08  which  corresponds  to  the 
position  in  the  profile  where  the  spanwise  oscillation  of  C  has  its 
maximum  amplitude.  It  is  also  evident  that  peaks  and  valleys  in 
C  are  lined  up  with  those  in  C ' . 

Further  downstream  at  x  =  9  cm  (figure  II. 8)  the  same  features 
are  evident  but  the  spanwise  spacing  between  peaks  and  valleys  has 
increased.  Using  the  spanwise  spacing  defined  by  the  central  part  of 

the  figure  gives  the  normalized  value  X/<S  =1.0  in  terms  of  the 

'  0) 

local  vorticity  thickness  (Konrad's  value).  Nearly  the  same 
dimensionless  value  is  obtained  at  the  upstream  station  (x  =  4  cm) 
when  an  average  value  for  X  is  calculated  from  the  three  wavelengths 
available  in  figure  II. 7.  Between  x  *  4  cm  and  x  *  9  cm  the 
change  in  scale  is  by  a  factor  of  1.9  in  5^  and,  evidently,  the 
same  for  the  spanwise  structure.  As  will  be  seen  in  figure  11.12, 
the  two  stations  are  just  upstream  and  downstream,  respectively,  of 
the  mixing  transition.  The  amplitude  of  the  spanwise  oscillations 
of  the  mean  concentration  is  not  strongly  affected  by  downstream 
position;  its  maximum  rras  variation  of  0.10  occurs  at  n  =  0.07. 

The  maximum  fluctuation  in  time,  whose  rms  value  is  c'  *  0.30  occurs 
at  the  same  value  of  n  but  as  noted  before,  has  little  spanwise 
variation. 

Some  insight  into  the  mean  spanwise  variation  can  be  obtained 
by  comparing  simultaneous  concentration  time  histories  from  the 
two  probes.  This  is  illustrated  in  figure  II. 9  where  the  concentration 


48 


histories  measured  by  the  two  probes  at  x  =  9  cm,  z  -  0  and 
z  =  0.64  are  shown  at  three  locations  across  the  layer.  At  the 
center  of  the  layer,  n  =  -  0.04,  large  concentration  fluctuations 
associated  with  the  large  scale  vortices  dominate  the  signal.  These 
large  fluctuations  appear  almost  simultaneously  in  time  indicating 
the  basic  two-dimensionality  of  the  large  scale  structures.  There  are, 
however,  differences  in  the  smaller  amplitude,  higher  frequency 
fluctuations  which  are  superimposed  on  the  large  ones.  At  the  edges, 

D  *  -  0.09  and  0.04,  there  is  in  general  a  correlation  between  large 
spikes  in  the  two  signals  but  their  shapes  and  deviations  are  not  so 
consistently  similar  as  at  the  center  of  the  layer  and  occasionally 
one  probe  shows  no  evidence  of  the  passage  of  the  vortex  seen  at 
the  other  probe  only  0.64  cm  away.  This  spacing  is  equal  to 
approximately  half  the  vorticity  thickness  at  this  location!  Clearly 
there  are  large  changes  in  the  internal  structure  of  the  large, 
spanwise  vortices.  As  will  also  be  confirmed  by  the  pictures  in 
Chapter  II. 3,  this  internal  structure  results  in  significant  distortions 
of  the  large  vortices  but  their  overall  spanwise  coherence  is  preserved. 

The  concentration  histories  were  also  used  to  calculate  the 
amount  of  chemical  reaction  product  (Konrad  1976).  The  mean  product 
formed  as  a  function  of  span  and  position  across  the  layer  is  plotted 
in  figure  11,10  for  x  *  4  and  9  cm.  Its  behavior  is  similar  to 
that  of  the  rms  of  the  concentration,  fluctuation,  with  a  180°  phase 
change  across  the  layer;  the  position  at  which  this  change  occurs  is 
r)  -  -  0.04,  different  from  the  value  found  for  the  rms. 


Wv».4^r. 


49 

The  product  thickness  defined  as 
00 

S~p  =  J  cp(y)  dy 

— oo 

is  plotted  versus  span  for  the  two  downstream  positions  in  figure  II. 11. 
The  low  amplitude  of  the  oscillations  results  from  the  observed  phase 
change  across  the  layer.  The  product  thickness  dependence  on  down¬ 
stream  distance  is  presented  in  figure  11,12  for  two  spanwise  positions 
z=0  and  z  =  0.64  cm.  Also  plotted  in  the  figure  are  the  measurements 
of  Konrad  (1976).  The  measurements  at  z  =  0  show  good  agreement 
with  those  of  Konrad,  in  particular,  the  mixing  transition  is 
observed  at  the  same  Reynolds  number.  The  measurements  at  z  =  0.64  cm 
show  a  smaller  amount  of  mixing. 

Comparison  was  also  made  between  the  result  of  Konrad  and  our  own 
with  the  time  exposure  pictures.  Shown  in  figure  11.13  is  a  time 
exposure  picture  together  with  a  plot  of  the  product  thickness  growth 
rate  6~  /(x  -  xfl)  versus  large  scale  Reynolds  number  Alfi^/v^ 
where  6^  is  given  by 

*  0.119  (x  -  x„) 

from  Konrad  (1976),  The  plot  and  picture  have  the  same  horizontal 
scale.  It  is  observed  that  the  mixing  transition  occurs  at  approxi¬ 
mately  the  downstream  location  where  some  of  the  streaks  in  the 
picture  merge,  resulting  in  the  scale  change. 
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II. 3  Constant  Density  Mixing  Layer 

11. 3.1  Apparatus  and  flow  visualization  technique 
The  water  mixing  layer  facility  designed  by  Breidenthal  (1978) 

was  used  in  flow  visualization  experiments  using  laser-induced 
fluorescence.  It  is  a  gravity  driven  apparatus  in  which  the  free 
stream  fluid  is  supplied  from  two  independent  reservoirs  located  one 
floor  above  the  test  section.  Although  the  facility  is  designed  to 
operate  at  higher  velocities  a  maximum  velocity  of  100  cm/s  was  used 
in  these  experiments . 

The  technique  of  laser-induced  fluorescence-,  which  was  developed 
in  this  laboratory  by  Professor  P.E.  Dimotakis,  was  used  for  flow 
visualization.  The  fluorescent  dye  fluorescein  was  dissolved  in  one 
free  stream.  A  plane  sheet  of  light  was  used  to  excite  its  fluorescence 
which  was  recorded  with  conventional  photographic  techniques,  The 
sheet  was  obtained  from  an  argon-ion  laser  and  had  a  thickness  less 
than  0.5  mm,  A  similar  system  has  been  used  by  Dewey  (1976). 

The  cross-sectional  flow  pictures  described  in  the  next  section 
were  obtained  with  the  plane  of  illumination  perpendicular  to  the 
free  stream  velocity.  Because  of  geometrical  limitations  in  the 
facility  the  photographic  and  motion  picture  camera  optical  axis 
formed  an  angle  of  approximately  30°  with  the  perpendicular  to  the 
visualized  plane,  which  resulted  in  some  distortion  of  the  pictures,  \ 

11. 3. 2  Flow  visualization  results  ; 

Two  cross-section  pictures  of  the  mixing  layer  at  a  Reynolds 

number  of  AU<$  /v  ■  3  x  103  are  shown  in  Figure  11,14.  The  flow  is 


labeled  fluid,  is  at  the  bottom.  Figure  11.14(a)  is  a  cross-section 
through  a  large  scale  primary  vortex  core  while  figure  11.14(b)  is  a 
cross-section  through  a  braid.  Even  though  the  technique  does  not 
visualize  vorticity,  the  presence  of  streamwise  counter-rotating 
vortices  is  inferred  from  the  swirled  interface. 

The  expected  changes  resulting  from  the  primary  vortex  structure 
are  readily  observed  in  these  pictures.  The  braid  is  a  region  where 
fresh  free  stream  fluids  meet,  resulting  in  a  sharp  contrast  change, 
while  in  the  core  there  is  a  region  of  mixed  fluid  separating  the 
free  streams.  The  streamwise  vortex  structure  also  shows  significant 
variations  in  the  two  regions.  In  a  cross-section  through  the  core, 
a  row  of  counter-rotating  vortices  stretching  across  the  span  is 
observed  at  both  sides  of  the  mixed  region  while  the  cross-section 
through  the  braid  has  a  single  row  of  streamwise  vortices.  (The 
properties  of  such  streamwise  vortices  in  a  braid  have  been  investigated 
by  Corcos,  1979.) 

The  peculiarities  of  the  streamwise  vortex  structure  and  its  inter¬ 
action  with  the  primary  vortices  was  further  investigated  with  a  high 
speed  motion  picture  obtained  at  a  Reynolds  number  AU6^/v  *  2.4  x  103. 
Three  sequences  from  this  movie,  shown  in  figures  11.15,  11.16  and  11.17 
will  be  discussed  in  the  following  paragraphs.  It  may  be  helpful  in 
following  the  discussion  to  refer  to  figures  11.18(a),  (b),  (c),  which 
were  constructed  on  the  basis  of  these  and  other  movie  sequences. 
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The  sequence  in  figure  11.15  shows  the  passage  of  a  primary 
vortex  through  the  visualized  plane.  In  the  first  two  pictures  the 
preceding  braid  is  seen  to  move  toward  the  low  speed  side,  and  the 
streamwise  vortices  follow  this  motion.  First  contact  of  the  primary 
vortex  core  with  the  plane  of  visualization  is  seen  in  the  third 
frame  (t  =  0.1  sec).  It  occurs  at  isolated  regions  along  the  span 
indicating  outflow  of  fluid  from  the  core  in  these  areas;  furthermore 
these  outflow  regions  correlate  well  with  the  counter-rotating  vortices 
at  the  low  speed  side,  in  particular  with  the  positions  where  their 
induced  velocity  field  is  directed  upwards. 

As  the  core  moves  through  the  visualized  plane  a  row  of  counter¬ 
rotating  vortices  appears  on  the  high  speed  side.  Vortex  pairs  on 
the  high  and  low  speed  sides  of  the  primary  vortex  are  aligned  with 
each  other  in  such  a  way  that  they  have  the  same  sign  of  vorticity  at 
the  same  spanwise  location.  Another  feature  of  the  streamwise  and 
primary  vortex  interaction  is  the  tendency  of  the  streamwise  vortices 
to  move  away  from  the  core.  This  trend  can  also  be  observed  in 
figure  11.14(a) . 

The  subsequent  evolution  shows  the  streamwise  vortices  at  the  low 
speed  side  disappearing  with  the  primary  vortex  core  leaving  behind 
a  region  of  outward  flow  from  the  core  (t  =  0.35)  i.e.,  normal  to  the 
picture  plane  upstream.  The  streamwise  vortices  at  the  high  speed 
side  continue  through  the  trailing  braid,  following  it  as  it  moves 
towards  the  low  speed  side. 
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The  outflow  region  observed  prior  to  the  appearance  of  the  stream- 
wise  vortices  on  the  high  speed  side  and  after  the  disappearance  of 
those  at  the  low  speed  side  indicates  that  the  streamwise  vortices 
bend  around  the  periphery  of  the  large  scale  vortex  core. 

No  evidence  was  found  of  the  continuation  of  the  streamwise 
vortices  along  the  periphery  of  the  primary  vortex  beyond  the  regions 
where  the  vortex  lines  are  perpendicular  to  the  plane  of  the  mixing 
layer,  which  suggests  that  the  vortices  in  each  counter-rctating  pair 
join  in  a  hairpin  configuration.  Although  direct  evidence  of  the 
presence  of  the  hairpins  was  not  obtained  this  conjecture  is 
consistent  with  the  observations  by  Breidenthal  (1978)  of  the  early 
development  of  three-dimensionality;  he  found  that  stretching  of  a 
primary  vortex  core  oscillation  (the  "wiggle")  leads  to  the  formation 
of  the  streamwise  vortices. 

A  slightly  different  evolution  is  observed  in  figure  11.16  where 
the  same  streamwise  vortices  at  the  low  speed  side  persist  through  two 
primary  vortices.  An  interesting  interaction  occurs  after  the  first 
primary  vortex  has  moved  through  the  visualized  plane  (t  =  0,4  to 
0.5  sec).  The  streamwise  vortices  at  the  low  speed  side  and  those  in 
the  intervening  braid  appear  to  retain  their  identity  until  the  second 
vortex  moves  through  the  plane.  This  hehavior  can  be  associated  with 
pairing  of  the  two  primary  vortices. 

The  sequence  in  figure  11,17  follows  the  streamwise  vortices 
between  neighboring  primary  vortices.  The  sequence  begins  at  a 
primary  vortex  core.  The  streamwise  vortices  are  first  seen  at 
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t  =  0.05,  in  particular,  a  counter-rotating  pair  on  the  high  speed 
side  slightly  to  the  right  of  the  center  of  the  picture.  As  the 
primary  vortex  moves  across  the  visualized  plane  the  two  vortices  in 
the  pair  separate.  The  lateral  motion  receeds  somewhat  as  the  vortices 
move  with  the  braid  to  the  low  speed  side.  The  process  accelerates  as 
the  streamwise  vortices  reach  the  low  speed  side  of  the  trailing 
primary  vortex.  This  behavior  is  associated  with  the  curvature 
imposed  on  the  path  of  the  streamwise  vortex  pair  by  the  primary 
vortex  velocity  field.  Furthermore,  this  curvature  induces  a 
precessional  motion  which  results  in  uneven  spacing  between  the 
streamwise  vortices  on  the  periphery  of  the  primary  vortex  and  there¬ 
fore  to  a  mutually  induced  radial  velocity  directed  outwards. 

These  induced  motions  outward  from  the  core  may  be  the  main 
contributions  to  the  observed  growth  between  pairings  of  the  primary 
vortices,  described  in  Part  I,  As  the  counter-rotating  pairs  move 
away  from  the  core,  their  induced  velocity  field  carries  mixed  fluid 
from  the  core  outward  while  it  brings  in  fresh  irrotational  fluid. 

Based  on  these  pictures  the  idealized  diagram  of  the  vortex  lines 
presented  in  figure  11,18  was  constructed.  A  streamwise  vortex 
segment  is  part  of  a  continuous  line  which  loops  back  and  forth 
between  adjacent  primary  vortices  as  shown  in  figure  11.18(a),  Thus 
on  successive  passes,  its  streamwise  direction  changes  sign,  forming 
the  vortex  pairs  that  are  seen  in  the  cross  section  views.  These 
warp-like  structures,  linking  adjacent  pairs  of  primary  vortices,  align 
themselves  with  each  other  so  that,  in  plan  view  they  appear  to  be 
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straight  streamwise  streaks.  Actually,  in  the  ideal  picture  the 
projections  onto  a  plan  view  would  not  really  be  continuous  but  would 
have  deflections,  above  and  below  the  primary  vortices  as  shown  in 
figure  11.18(c),  due  to  the  induced  precessional  motions  described 
above.  On  shadow  and  schlieren  pictures,  these  details,  superimposed 
on  other  features  from  the  internal  structure  of  the  primary  vortex, 
are  obscured. 

In  figure  11.19  is  a  sequence  from  a  high  speed  motion  picture 
obtained  at  a  Reynolds  number  AUS^/v  =  104  which,  according  to  the 
results  of  Breidenthal  (1978),  is  downstream  of  the  mixing  transition. 
The  reasons  for  the  poorer  resolution  in  this  sequence  are  not  entirely 
clear,  but  it  is  probably  simply  due  to  poorer  focussing  resulting 
from  the  wider  aperture  that  had  to  be  used  at  this  higher  flow  speed. 
The  figure  shows  the  passage  of  a  braid  through  the  visualized  plane. 
The  streamwise  vortices  can  be  followed  along  the  braid.  Their 
spacing  Is  smaller  than  in  previous  pictures  since  the  higher  velocity 
results  in  a  thinner  initial  shear  layer,  which  determines  their 
spacing. 
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II, A  Discussion  and  Conclusions 

Perhaps  the  most  striking  result  of  this  investigation  is  the 
elucidation  of  the  streamwise  vortex  structure  and  its  spatial 
relation  and  interaction  with  the  primary  spanwise  vortices.  An 
idealized  sketch  is  given  in  figure  11,18.  Unlike  previous  interpre¬ 
tations  of  the  streak  patterns,  that  each  streamwise  vortex  is 
unidirectional,  it  may  now  be  seen  that  it  is  part  of  a  warped  vortex 
line  threading  its  way  up-  and  downstream  between  any  two  adjacent 
primary,  spanwise  vortices,  thus  changing  its  streamwise  directional 
sign  on  each  pass.  This  structure  is  first  established  after  the 
Kelvin-Helmholtz  instability  has  developed  into  spanwise  coherent 
vortices,  possibly  as  early  as  the  first  pairing.  The  ’’wiggle" 
described  by  Breidenthal  (1978)  is  its  first  manifestation. 

As  defined  by  the  time  exposure  pictures,  the  mean  onset  Reynolds 
number  was  found  to  decrease  as  the  velocity  ratio  increase,  i.e.,  as 
the  mixing  layer  approaches  a  wake.  These  results  indicate  that  an 
initial  profile  containing  a  wake  component  is  more  unstable  to  three 
dimensional  disturbances.  Also  determined  from  the  time  exposures,  the 
mean  spanwise  spacing  of  the  streaks,  normalized  by  the  mean  local 
spacing  of  the  primary  vortices,  is  independent  of  velocity  ratio  and 
equal  to  0,67,  This  is  equal  to  the  most  amplified  spanwise  wave 
length  found  by  Pierrehumbert  (1980)  in  his  computation  of  the 
translative  instability  of  a  row  of  Stuart  vortices.  The  scatter  in  the 
measurements  is  consistent  with  the  broad  amplification  curve. 

As  pointed  out  by  Pierrehumbert  the 
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shallow  peak  is  also  consistent  with  the  fact  that,  once  formed,  the 
structure  preserves  its  spanwise  scale  for  some  distance  downstream, 
through  one  or  two  additional  pairings. 

The  proposed  warped  vortex  line  structure  was  inferred  mainly 
from  a  frame-by-frame  study  of  cross  sectional  movies.  It  is  also 
consistent  with  the  concentration  measurements,  which  show  that, 
associated  with  this  structure,  there  is  a  well  defined  spanwise 
entrainment  pattern.  The  fluid  from  the  high  and  low  speed  side  is 
entrained  at  different  spanwise  positions.  The  fluids  are  mixed  in 
the  primary  vortex  core  and  the  resulting  mixture  is  pulled  outwards 
at  different  spanwise  positions.  The  rms  of  the  concentration 
fluctuation  and  the  mean  mixing  product  show  similar  features 
consistent  with  the  pattern  described  above. 

The  observed  pattern  should  also  have  an  effect  on  the  spanwise 
velocity  distribution.  In  regions  where  high  speed  fluid  is  entrained 
a  larger  momentum  is  expected.  Measurements  by  Jimenez  et  al.  (1979) 
show  precisely  this  effect,  although  their  interpretation  is  somewhat 
different  from  ours. 

The  concentration  measurements  also  show  that  eventually  the 
scale  of  the  secondary  structure  increases.  The  mean  position  for  this 
has  not  been  precisely  determined  but  both  the  concentration  measure¬ 
ments  and  the  time  exposures  indicate  it  to  be  near  the  end  of  the 
mixing  transition.  The  streaky  pattern  is  still  seen  on  instantaneous 
pictures  far  downstream  of  the  mixing  transition  (e.g.,  Konrad,  1976) 
and  on  high  speed  schlieren  motion  pictures  which  we  obtained. 
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Downstream  of  the  transition  the  pattern  does  not  come  through  on  time 
exposures  because  it  varies  from  frame  to  frame,  apparently  with  the 
scale  of  the  primary  vortices , which  now  have  a  broad  distribution. 

The  mean  value  of  the  spanwise  spacing  now  scales  with  the  local 
vorticity  thickness,  specifically  \/5^  =  1-0.  (Jimenez  et  al. 
proposed  A/6  ~  1  to  1.25  on  the  basis  of  their  spanwise  measurements 

U) 

of  the  velocity  field.)  It  is  also  interesting  to  compare  this  result 

with  velocity  correlation  measurements  at  high  Reynolds  number  by 

Browand  and  Troutt  (1980)  who  found  a  large  value  of  correlation 

coefficient  (  >  80%)  for  spanwise  spacing  smaller  than  6  ,  dropping 

ui 

to  40%  for  2.36^.  Degradation  of  the  correlation  is  possibly  related 
to  the  secondary  vortex  structure. 

Finally  we  note  an  interesting  connection  with  the  observations, 
in  Part  I,  of  the  growth  of  the  primary  vortices  during  their  life¬ 
times,  as  seen  in  edge  view.  Much  of  that  growth  may  be  the  result  of 
the  induced  motions,  away  from  the  primary  vortex  cores,  of  the 
warped,  secondary  vortices. 

Whether  there  is  need  for  additional  instabilities  to  account  for 
the  mixing  transition  is  not  clear.  The  secondary  vortices  themselves 
increase  the  interface  considerably,  as  seen  in  the  cross-sectional 
pictures.  Furthermore,  after  the  vorticity  is  increased  by  stretching 
along  the  braids  and  periphery  of  the  cores,  amalgamation  of  the 
primary  vortices  results  in  accumulation  of  that  stretched  vorticity 
into  the  newly  formed  core.  Consequently,  entrained  but  unmixed  fluid 
in  the  spanwise  vortex  cores  in  the  initial,  two-dimensional  region 
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will  be  mixed  progressively  after  the  appearance  of  the  secondary 
vortex  structure  and  a  few  amalgamations.  Comparison  of  the  mixing 
transition  Reynolds  number  with  that  of  the  onset  of  the  secondary 
structure  indicates  that  (for  s  *  7  and  r  =  0.38)  it 
encompasses,  on  average,  2.5  amalgamations  by  pairing. 


Epilogue 


The  coherent  structure  of  turbulent  mixing  layers  has  been 
investigated.  The  results  of  several  experiments  and  theoretical 
developments  dealing  with  the  primary  spanwise  and  the  secondary 
streamwise  vortex  structures  have  been  presented.  It  is  hoped  that 
these  provide  some  further  insights  into  the  physical  processes 
occurring  in  a  turbulent  mixing  layer  and  will  motivate  additional 
research  in  this  challenging  field  of  Fluid  Mechanics. 
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xt/0i 

Bradshaw  (1966) 

1000 

Breidenthal  (1978) 

200 

800 

349 

250 

35 

Hussain  and  Zedan 

204 

310 

110 

(1978) 

185 

600-400 

214-150 

Large  Scale 

Self -similarity 

_ _ _ _  -  -  —  . . -  .  J 

180 

560 

200 

Table  1.1  Summary  transition  data  for  constant  density 
mixing  layers  with  initial  laminar  boundary 
layers. 


n 

0 

L/X 

Tearing 

3 

.1617 

.7975 

1.1983 

.4017 

Pairing 

2 

.2760 

.6602 

1.3301 

.6731 

Tripling 

3/2 

.4363 

.4840 

1.4947 

1.0216 

Table  1.2  Summary  of  theoretical  results. 
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Figure  1.1  Mixing  layer  parameters  and  coordinate  system. 


Figure  1.2  Shadowgraph  optical  system. 
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Figure  1.6  Mean  circulation  as  a  function  of  downstream  position. 
{ - ),  least  squared  fit  to  measured  points  with  x  >  4  cm. 


Figure  1.7  Normalized  rms  deviation  of  the  circulation  as  a  function 
of  downstream  position.  ( - ),  average  value  x  >  4  cm. 


Figure  1.8  Vortex  circulation  probability  distribution.  Histogram 
present  measurements;  (---),  Brown  and  Roshko  (1974);  (-•-),  Takaki 
and  Kovasznav  (1978);  ( - ),  Theory  chapter  1.3. 


Figure  1.12  Eulerian  description  (x  fixed):  (a)  Vortex  circulation 
probability  distribution  p(£)  (Eq.  1.3.14);  (b)  Vortex  formation 
probability  distribution  pp(C)  (Eq.  1.3.16);  (c)  Vortex  annihilation 

probability  distribution  PACS)  (Eq.  1.3.16). 


Figure  1.13  Lagrangian  description  (£  fixed):  (a)  Vortex  position 
probability  distribution  p-CO  (Eq.  1.3.19);  (b)  Vortex  formation 
probability  distribution  Pp(G)  (Eq.  1.3.20);  (c)  Vortex  annihilation 
probability  distribution  p^(£)  (Eq.  1.3.20). 


Photographic 

Plate 


Figure  II. 1  Schlieren  optical  system. 
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Figure  II. 2  Schlieren  pictures  of  the  helium-nitrogen  mixing  Inver.  ^ 

Side  view  upper  half,  plan  view  lower  half.  Uj  =  496  cm/s,  r  =  0.5,  j 

s=7,p0=4atm.  3 

(a)  1  ps  exposure.  j 

(b)  2  second  exposure. 
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(a) 

Figure  II. 3  Same  as  figure  II. 2  at  conditions  Uj  =  1470  cm/s, 
r  =  0.3,  s  =  7,  p0  =  4  atm, 

(a)  1  ys  exposure. 

(b)  2  second  exposure. 
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Figure  II. 5  Normalized  streak  spacing  as  a  function  of  AU/Uc,  s  » 


igure  II. 6  Strearawise  evolution  of  mean  concentration  profiles  at  two 
panwise  positions:  ©  ,  z  -  0;  □  ,  z  -  0.64  cm,  (Measured  points  in  the 
rofile  are  connected  by  straight  line  segments.) 
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(b) 

Figure  II. 7  Spanwise  variation  of  the  mean  concentration  (a)  and  Its 
rms  fluctuation  (b)  at  various  locations  in  the  profile  at  x  -  4  cm. 


variation  of  mean  product. 
4  cm. 
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Figure  11.14  Cross-section  pictures  of  the  mixing  layer  in  water. 
Re<t  *  3  *  103. 

(a)  Cross-section  through  the  core. 

(b)  Cross-section  through  the  braid. 
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Figure  11.16  Cross-section  movie  sequence  of  a  primary  vortex 
amalgamation  Re^  ■  2,4  x  10 3 , 


